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We rediscuss the question of what fraction of the population of MACHOs in the 
halo of the Milky Way might be in the form of halo dwarfs, because recent deter- 
minations indicate masses of the MACHOs, which are typical for such stars. For 
this purpose we have analyzed the Third Catalogue of Nearby Stars and identified 
the subdwarfs and a high velocity white dwarf in the solar neighbourhood. The 
local mass density of these stars is 1.5-10 - 4 A4q pc -3 , which is only 3% of the 
current estimate of the local mass density of the MACHO population. We com- 
pare the local density of subdwarfs with constraints set by HST observations of 
distant red dwarfs. Using models of the stellar halo with density laws that fall off 
like r~ a , a = 3.5 to 4, we find that the HST constraints can only be matched, if 
we assume that the stellar halo is flattened with an axial ratio of about 0.6. The 
non— detection of the analogs of MACHOs in the solar neighbourhood allows to set 
an upper limit to the luminosity of MACHOs of Mg > 21 magnitudes. 

1 Introduction 

The MACHO § and EROsi collaborations have reported new results of their 
campaigns of observing micro-lensing events towards the LMC. These indicate 
that MACHOS have masses typically of half a solar mass and contribute about 
one half to the mass budget of the halo of the Milky Way (see the article by 
Bennett in this volume). n Such masses are typical for red and white dwarfs. 
Bahcall and collaboratorsHQ have made by very deep star counts based on HST 
data in situ measurements of the space density of such stars in regions, where 
the micro-lenses are expected to be physically located. Their conclusion is that 
halo dwarfs contribute only a few percent to the mass of the halo and are thus 
unlikely MACHO candidates. Graff & FreeseB (see also Freese's article in this 
volume) infer from the same data that the local mass density of halo dwarfs 
is even less than one percent of the combined local densities of the dark and 
stellar halos. 

On the other hand, halo dwarfs are directly observed in the solar neigh- 
bourhood. LiebertB discusses the results of the USNO parallax programme 
for red dwarfs and finds that the number of halo dwarfs in that sample is 
consistent with the HST data. We have followed a complementary approach. 
Recently the Third Catalog of Nearby Stars (CNS3) has been completed at 
the Astronomisches Rechen-Institut. This provides the now most complete 
inventory of the solar neighbourhood up to a distance of 25 pc from the Sun 
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Figure 1: Colour-Magnitude diagram of nearby stars with um v < 0.3 mag. 

and allows a new determination of the local density of halo stars. 
2 Subdwarfs in the CNS3 

We have searched the CNS3 for halo dwarfs. These reveal themselves by their 
position on the subdwarf sequence in the colour-magnitude diagram (CMD), 
their low metallicity, and their high space velocities. Unfortunately, there is in 
this metallicity range a broad overlap of the comparatively few halo stars and 
the much more abundant thick disk stars. We have therefore adopted a very 
conservative search strategy in order to isolate a true halo population, even if 
this will lead to an undersampling of the halo population. We have examined 
first the stars on the subdwarf sequence, which can be clearly seen in the CMD 
shown in Fig. 1. The reliability of the photometry and the parallaxes of each 
star has been checked carefully, which led to the omission of quite a number 



Table 1: Subdwarfs in the CNS3. 
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of stars with inaccurate parallaxes. In this way we identified G1191, G1781, 
WO 9722 = LHS 64, GJ1062, LHS 375, and LHS 3409 as likely halo stars. A 
number of stars in the CNS3 are classified as subdwarfs by their spectral type. 
This enables us to search for subdwarfs in that part of the CMD, V — I < 
1.5, where the subdwarf sequence comes very close to the main sequence. We 
include GJ 1064A+B and G153A+B, which have reliable parallaxes, into our 
list on the basis of this criterion. All stars have decidedly large space velocities 
(cf. Table 1 and Fig. 2). In addition to these stars we found 4 stars, G1158, 
G1299, G1451A, and the white dwarf Gl 699.1 (spectral type DA7), which do 
not lie very prominently on the subdwarf sequence, but have space velocities, 
which clearly identify them as halo stars. 

All the stars in our list have very low metallicities. Leggetti has deter- 
mined the metallicities of a large number of nearby M dwarfs from their posi- 
tion in the (J-H)-(H-K) two infrared-colour diagram. According to that deter- 
mination G1191, G1781, GJ1062, and G1299 have metallicities [Fe/H] < -1. 
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LHS 375 has been classified as a cold subdwarf by Ruiz and Anguita □. The 
metallicities of G153A, G115JS, G1451, and_GJ1064A, B have been taken from 
the compilation by Taylor EE Reid et al. B have classified GJ 1062, G1781B, 
WO 9222, and LHS 3409 as likely subdwarfs on the basis of their spectra. Jones 
et al.E3 confirm the low metallicity of G1299. There are many apparently low 
metalliciiy stars in the catalogue, which must be thick disk stars. For instance 
LeggettQ has assigned to 20 stars in the CNS3 a metallicity of [Fe/Hj < — 1, 
of which only 4 have been finally included in our sample. 

In Table 1 we give the list of stars selected out of the CNS3 as described 
above, of which we are now reasonably certain that they are genuine halo stars. 
The velocity components of the subdwarfs given in columns 6 to 8 in Table 1 are 
referred to the Sun. In order to reduce them to the LSR the solar motion (U, V, 
W) Q = (+9, +12, +7) km/s has to be added. Note that U points towards the 
galactic center. The velocity distribution of our sample is shown as a Toomre 
diagram in Fig. 2. It is consistent with a distribution centered on V = -220 
km/s. No clustering near V = 40 km/s, which would indicate a contamination 
by thick disk stars, is detected. The dispersion of \/U 2 + W 2 is 150 ± 3D km/s, 
which is within statistical errors consistent with other determinationsl!3 of the 
kinematical parameters of halo stars. The masses given in the last column of 
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Figure 2: Toomre diagram of the nearby halo dwarfs. The dashed lines indicate lines of 
constant kinetic energy of the halo stars. 

Table 1 have been assigned using the mass-to-luminosity relation-calculated 
theoretically for a metallicity of [Fe/H] = -1.35 by Alexander et allij. 



a G1781 is a variable star and shows H a emission. However, it is a spectroscopic binary, 
which, in our view, accounts for this phenomenon. 
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3 Results and Discussion 



3.1 Local density 

The local mass density of halo stars can be estimated in various ways. In Fig. 3 
we show density estimates, which have been calculated from the cumulative 
distribution of the stars in our sample, i.e. by adding up the masses of stars 
within given distances and dividing by the corresponding spherical volumes. 
The stars within 10 pc give rather high density estimates of the order of 1CP 3 
M.q pc~ 3 . Even though Wielen & JahreifiEj found similar values, when eval- 
uating the second edition of the Gliese catalogue, they seem unrealistic to us. 
We find a plateau in the density run, when we consider the stars within the 20 
pc sphere (cf. Fig. 3). Beyond that the density drops off, probably because the 
catalogue becomes incomplete. The most likely value of the local mass density 
lies according to this determination in the range 1.5 to 110 -4 A4q pc~ 3 . The 
local density of the dark haloEEl is about 910 -3 M.& pc -3 (see also Gates' 
article in this volume), so that the local density of halo stars is about 1.7% 
of the halo density or 3% of the local MACHO density as determined by the 
MACHO collaboration. We have broken down our sample of subdwarfs with 
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Figure 3: Mass density of nearby halo dwarfs as function of the distance from the Sun. 
Within 10 pc only the positions of the stars are indicated. The errorbars indicate the 

statistical uncerntaintics. 



respect to absolute magnitudes. The resulting 'luminosity function' has been 
compared with the subdwarf luminosity function obtained by Dahn et al. E3 
for stars with My > 10 mag. Both luminosity functions agree well within 
statistical errors. 



5 



3.2 Matching the constraints by HST data. 

It is interesting to confront the density of halo stars derived here with density 
estimates based on the HST data. Graff & Freesetl find in a detailed analysis 
of earlier HST data la a value about five times less than our estimate, which 
is partly due to their restriction to red dwarfs redder than V — I ~ 2 mag. 
Whereas these authors tried to infer the local density of halo stars from the star 
counts, we reverse the approach and make predictions of the expected number 
of stars in the star count field. Very recently Flynn et al. □ have reported an 
analysis of the Hubble Deep Field (HDF). The HDF has a size of 4.4 square 
arcminutes and is located at I = 126°, b — 55°. The limiting magnitudes are 
18 to 26.3 mag in the /-band. 

In the following we assume that the population of halo stars follows a 
density law of the form 



-=-4^* ) m 

with an arbitrarily chosen core radius of r c = 1 kpc. Such density laws-Are 
typical for other tracerSj-af the halo population, such as RR Lyr stars ta or 
horizontal branch stars tL3, in the Milky Way or external galaxies E3. The 
exponent a lies typically in the range 3.5 to 4. Note that the dark halo is 
usually described by a density law of the form as in equation (1) with a = 
2. If the subdwarfs are observed locally at a density Vh Q , one predicts by 
integrating along the line of sight 

N = ^ j(T 52 ( rl+rlJ-t^lcosb )^ 3 (2) 

stars in the star count field. Q is the angular area of the field and the minimum 
and maximum distances are determined by the limiting magnitudes. We con- 
centrate first on the colour range V — I = 1.8 to 3.5 mag. We have 4 stars in 
that range within 16.5 pc, where we believe our sample to be reasonably com- 
plete. If we project these into the cone towards the HDF, we obtain the star 
numbers summarized in Table 2. The limiting distances have been calculated 
for each star individually, 



A mm,max 



1 Q0.2(J mta , max +(V r -I)-Mv+5) ^ 

where, in order to avoid confusion with disk stars in the HDF, we adopt a lower 
limit of I m in = 24.6 mag. As can be seen from Table 2 we predict 7 to 11 stars 
in the HDF, whereas Flynn et al. □ have actually dectected no star. Despite 
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Tabic 2: Predicted Number of Stars in the HDF and the GS. 



V -I > 1.8 
a c/a hjjdf 


V -I < 1.8 
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4 0.6 2 
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2 32 



the low number of stars involved, this discrepancy seems to be statistically 
significant. Several explanations might account for this discrepancy. First, the 
stellar halo might be much more irregular and lump}? thaji .previously thought. 
Second, the stellar halo is almost certainly flattened Il§t2l23 with an axial ratio 



around c/a w 0.6. If we take such a flattening into account in the halo model, 

( r c+ R % Y 

Vh = Vh ®\ 1 * + R* + z y(c/a)*) ' (4) 

where R, z denote cylindric coordinates, we predict star numbers in the HDF, 
which are statistically consistent with no star seen by Flynn et alu. 

The star counts in the HDF in the colour range V — I < 1.8 mag can be 
interpreted in a similar way. Disk stars in this colour range are so bright that 
they would have to lie several kpc above the midplane to appear fainter than 
I ~ 22 mag. Thus in order to avoid confusion with disk stars in the HDF, we 
consider only stars fainter than I = 22 mag. In our sample there are 4 stars 
within 20 pc in this colour range. Their V — I colours cluster around 0.9 and 
the white dwarf has V — I = 0.5. Thus we define a colour range 0.5 < V — I < 
1.8, in which we compare the predictions with actually observed numbers of 
stars. If the stars of our sample are projected into the cone towards the HDF 
using spherical halo models, we predict 7 to 17 stars depending on the model, 
whereas Flynn at al.Q have observed 4 stars. If we assume again a flattening 
of c/a = 0.6, the number of expected stars is consistent with the observed 
number. 

Gould et al.cd are presently investigating another field observed by HST. 
This is the so called Groth strip (GS), which has an angular size of 114 square 
arcminutes and is located at-i = 96° and b = 60°. Its limiting magnitude is 
/ = 23.9 mag. Gould et al.cd have very kindly made their preliminary data 
available to us, so that we were able to make star counts in this field. Due to 
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the limiting magnitude of I — 23.9 mag even the faintest stars in the GS with 
colours redder than V — I = 1.8 could be disk stars. Thus we focus on the 
colour range V — I < 1.8 and adopt the same halo zone as in the HDF. The 
GS has 66 stars in this zone. If we project the local subdwarfs into the cone 
towards the GS, we obtain the numbers of predicted stars summarized in the 
last column of Table 2. Again a flattened halo model is a better fit to the data. 



3.3 Constraints on the absolute magnitude of MACHOs. 

Like previous authors we have not detected the analogs to MACHOs in the solar 
neighbourhood. This raises the question why. All subdwarfs in our sample 
appear in the LHS high-proper motion star catalogue. Being halo objects, 
MACHOs in the solar neighbourhood must have similar proper motions. The 
fact that they have not been found in high-proper motioruaurveys allows to 
set an upper limit to their luminosity. The LHS catalogueE 2 ], for instance, is 
claimed to be complete for stars down to apparent magnitude B — 21 mag and 
with proper motions 2. 5/a > /i > 0. 5/a. For stars brighter than B = 10 mag 
it is claimed that all stars with proper motions fj, > 0. 3/a are contained in 
the catalogue. We assume that the MACHOs are homogeneously distributed 
around the Sun and have a gaussian velocity distribution, 

dn = to W2 3 ex P -^(U 2 + (V- Vf + W 2 )d 3 vd 3 r , (5) 

centered on V = -220 km/s and with a velocity dispersion typical for halo 
objects, ctm = V/y/2. Integrating now over all radial velocities and the proper 
motions according to the specifications of the LHS we can determine the radius 
r2 of the sphere out of which one would expect say 2 MACHOs in the LHS, 
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where Iq denotes the Bessel function. According to Poisson statistics two 
expected MACHOs would be still consistent with no MACHO actually seen 
in the LHS survey. In Table 3, we give results of numerical integrations of 
equation (6) assuming a local density of MACHOs of vm q = 0.5 • 0.01 A4q 
pc~ 3 / 0.5 Mq as determined by the MACHO collaboration. A lower limit 
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of the absolute magnitude of the MACHOs is then given by Mb — 21 — 
51ogr2 + 5, because MACHOs more distant than will not show up in the 
LHS, since they are fainter than the apparent magnitude limit of the LHS, 
B = 21 mag. From Table 3 we conclude that MACHOs are fainter than Mb 
— 21.2 mag. The southern hemisphere (5 < —30°) and the galactic belt are 
not sampled by the LHS survey. Thus one could argue that only the region 
b > 30° is properly sampled. This would reduce the numbers in the second 
column of Table 3 by a factor of 0.25 and shift the lower limit of the absolute 
magnitude of the MACHOs to M B = 20.6 mag. If the distribution of MACHOs 
is irregular and lumpy, the micro-lensing results could mimic a too high mean 
density of MACHOs. If we assume a local MACHO density of only 10% of the 
value determined by the MACHO collaboration, the lower limit of the absolute 
magnitude of MACHOs would be Mb = 20.1 mag. 



Table 3: Estimated absolute magnitudes of MACHOs. 
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If the MACHOs were identified with faint red dwarfs, their masses would 
be less than 0.1 A4q in contradiction to the estimate by the MACHO collab- 
oration. As an alternative very faint white dwarfs have been discussed in the 
literature as MACHO candidates^. If we use the theoretically calculated cool- 
ing sequence of a DA white dwarf by WoodEj and extrapolate B — V to 2 mag, 
we estimate cooling times of 11 to 13 Gyrs for white dwarfs as faint as Mb = 
20.1 to 21.2 mag. This would make them as old as globular clusters and raises 
the question, whether this allows for reasonable lifetimes of the progenitors of 
the white dwarfs. 
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